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, We apply the potential force estimation method to artificial time series of market price produced 

C ' by a deterministic dealer model. We find that dealers' feedback of linear prediction of market price 

based on the latest mean price changes plays the central role in the market's potential force. When 
markets are dominated by dealers with positive feedback the resulting potential force is repulsive, 

' while the effect of negative feedback enhances the attractive potential force. 
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1. INTRODUCTION 



Improvement of time resolution of market data owing to the development of computer 
technology has been led to many new discoveries in econophysics. Empirical laws have been 
established in foreign exchange markets and stock markets such as the power laws in price 



change distributions, the long-tai 
interval distributions, and so on 



s in volatility auto-correlation, non-Poissonian transaction 



H [21 01 • One of the goals of econophysics study may be to 
clarify the origin of these empirical laws. As the market prices are actually determined by 
the traders' transactions, it is natural to expect the origin of these empirical laws to be the 
dealers' actions. In our previous study of deterministic agent model which we call the dealer 
model, we successfully reproduced most of the basic empirical laws j4|[5|Q]. 

Recently, M.Takayasu and her coworkers introduced a new time series analysis method 
to observe non-trivial time dependent potential forces in market prices It is shown that 
this potential force can be observed universally in financial markets, and the most of basic 
empirical laws are also satisfied in this formulation. In this paper we numerically observe 
potential forces of the dealer model and clarify an origin of market potentials. 



2. POTENTIAL FORCES OBSERVED IN FOREIGN EXCHANGE MARKETS 

In this section we briefly show an example of results of the potential analysis to a foreign 
exchange market data for review of the method and for comparison with the dealer model 
in the next section (Fig{T]). The figure located in the center is the optimal moving average 
price P{u) defined by the tick time u updated by every transaction. The optimal moving 
average is obtained by removing uncorrelated noise component 9| from dollar- yen bid price 
during about one week from January 12th 1999 to January 14th 1999. The intervals marked 
by A, B, C and D consist of 2000 ticks and we estimate potential functions directly from 
these ticks following the method^]. And in the sub-figures A-1, B-1, C-1 and D-1, we plot 
{P{u) — Pm{u)) vs (M — l)f/(P) for each intervals. Here, the super- moving average Pm is 

defined as Pm = P{u — k) where the size parameter M is 16. U (P) denotes the potential 

fc=0 

function, and these variables are assumed to satisfy the following equation with a stochastic 
noise term, Fiu): 
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P(n+ 1) = P{u) 



d_ 
dP 



U{P) \p=Piu)-PMiu) +F{u) 



(1) 



U{P) 




(2) 



As known from these figures we can approximate the plots by quadric potentials and 
estimate the curvature of the potential, b{u), as a function of time. 

In the intervals A and B the curvature b{u) is positive and the price fluctuation is attracted 
by the center of price Pm- The diffusion of price fluctuation is slower than a simple random 
walk case with the same volatility. On the other hand in the intervals C and D, the potential 
coefficient b{u) is negative and the diffusion is faster. As shown the figures, we can clear 
time dependent potentials in market prices. It is confirmed that for pure random walk data 
or randomly shuffled data the estimated potential coefficient distributes around with the 
standard deviation about 0.2. 

3. COMPARISON OF POTENTIALS WITH THE DETERMINISTIC DEALER 
MODEL 

In this section we apply the potential analysis to the time series produced by the de- 
terministic dealer model. In this model we assume that this artificial market consists of 

= 300 dealers and every dealer changes his bid price Pi{s) at time s described by the 
following equation. 



The parameter q is the amount of change at each time step: if i-th dealer is a buyer, he 
monotonically raises his price until he can buy while if he is a seller he reduces his price 
until he can sell. We set Cj uniform random numbers given initially to represent dealers' 
character, namely, a dealer who has big Cj is regarded as quick-tempered. The second term 




(3) 



Here, ctj represents the i-th dealer's position. 




+1 buyer 



1 seller 
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di < AP > M represents the foreseeing effect defined in the following. 
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M 



M(M + 1) 



^(M-(fc-l))AP(n-(fc-l)), 



(4) 



fc=i 



where AP{u — k) = P{u — k) — P{u — {k + 1)) is price change before k ticks. This is a 
weighted moving average of price changes for past M ticks with bigger weights to newer 
price changes. For simplicity we set the value of di constant for all dealers, that is, di = d. 
When d is positive, the dealers are trend followers who predict the near future market price 
change to be proportional to the latest price change. When d is negative, the dealers are 
contrarians who predict in the opposite way, for example, when the market price is going up 
he predicts that the price will go down in the near future. It is possible to tune the values 
of parameters so that some of the empirical laws of financial markets are reproduced by this 
formulation [slj^. 

In this model the j-th dealer's ask price is given by Pj{s) + L assuming that the spread is 
a constant, independent of the dealers. In this artificial market a transaction occurs when 
a bid price and an ask price match. This condition is described as follows: 



where max{pj(s)} is the maximal value of bid prices and min{pj(s) + L} is the minimal 
value of ask prices. Transaction occurs between the pair of dealers who quote the highest 
bid price and lowest ask price, and the market price P{u) is defined by the middle value, 
[max{pj} + min{pi + L}]/2. The dealers are assumed to trade unit volume and after the trade 
the dealers change their positions, so the buyer becomes a seller and the seller becomes a 
buyer. In other words, the sign of at is reversed when the i-th dealer is involved in a 
transaction. By repeating these deterministic evolution rules this artificial market keeps 
producing apparently random market prices continuously. 

We apply the potential analysis method previously described to the resulting artificial 
market price data. For any value of the foreseeing parameter d we can observe non-trivial 
potential forces and the value of potential coefficient can be estimated. Here, we show four 
typical cases marked A', B', C, D' in Fig|2l In these simulations all parameters except d 
are identical. In the interval of Figl2] A', the parameter d is negative, that is, dealers are 
contrarians and we find a stable potential as shown in Fig|2]A'-l. When the parameter d 
equals zero, the observed potential function is nearly fiat as shown in Figj2]C'-l. When the 



max{pj(s)} < mm{pj{s) + L}, 



(5) 
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parameter d is greater than zero, the dealers are trend-followers and we observe an unstable 
potential as shown in Fig|2]C'-l and D'-l. From these results we understand that the trend- 
following effect directly contributes to the unstable potential and the contrarian effect is 
responsible for the stable potential. 

In FigJ3]we plot the estimated value of the potential coefficient b* taken average over b{u) 
for 100,000 ticks as a function of the foreseeing parameter d. All the points are clearly on a 
straight line b* = 0.2 — 0.86d. When d = 0, we can observe very weak attractive potential 
force, b* = 0.2 by taking average over long numerical simulation. We conjecture that this 
attractive force is not a numerical error but it is caused by the transaction rule. In our 
model a buyer becomes a seller and a seller becomes a buyer after the transaction. Due to 
this effect of changing roles between two dealers, it pulls back the market price to the mean 
values of the past prices. It can be described by the random walk in the attractive force. 

4. CONCLUSION AND DISCUSSION 

In this paper we showed by numerical simulations that dealers' foreseeing effect is the 
source of the market potential force observed in real financial markets. As a result, we 
have found that the coefficient of the potential function and the dealers' foreseeing action is 
deeply related. When the dealers are trend-followers the potential coefficient is negative and 
the market is unstable, and when the dealers are contrarians the coefficient is positive and 
the market is stable. In real markets the coefficient value is always fluctuating as typically 
shown in Fig.l. This implies that the prediction strategy of dealers in real markets may not 
be fixed to either trend-followers or contrarians, but they change the strategy from time to 
time. 

For real market dealers it is said that to know other dealers' strategy or the whole market 
atmosphere is very important. It may become possible to tell by using the coefficient b 
whether the market is dominated by trend-followers or contrarians in real time. 



5 



Acknowledgments 



[1] T. Mizuno, S. Kurihara, M. Takayasu and H.Takayasu, Physica A 324 (2003) 296-302. 
[2] P. Gopikrishnan, M.Meyer, L. A. N. Amaral, V. Plerou, and H. E. Stanley, Eur. Phys. J. B 3 
(1998) 139-140. 

[3] M. Takayasu, H.Takayasu, M.P.Okazaki, Proceedings of "Empirical Science of Financial Fluc- 
tuations" in Tokyo, (edited by H. Takayasu, Springer, 18-26, 2001) 

[4] H.Takayasu, H.Miura, H.Hirabayashi and K.Hamada, Physica A 184 (1992) 127-134. 

[5] A.Sato and H.Takayasu, Physica A 250 (1998) 231-252. 

[6] K.Yamada H.Takayasu and M.Takayasu, Physica A 382 (2007) 340-346. 

[7] M.Takayasu, T.Mizuno and H.Takayasu, Physica A, 370 (2006) 91-97. 

[8] L. Bacheher, Ann. Sci. Ec. Norm. Super, (iii) 17 (1990) 21 (Translation: Cootner, 1964) 

[9] Takaaki Ohnishi, Takayuki Mizuno, Kazuyuki Aihara, Misako Takayasu and Hideki Takayasu 
Proceedings of gPractical Fruits of Econophysicsh (edited by H. Takayasu, Springer, 62-66, 2005) 



6 



X 



114 — 



0) 

u 

Cl 112 



110 




X -2 



_L 




5 1 O 

Time [ticks] 



X 



_L 



1 5x1 O 




P-Pm 



FIG. 1: Calculated potentials by Takayasu's method in a real market. The centered time series 
is noise free prices so-called optimal moving average during about two days from January 12th 
1999 to January 14th 1999. The ranges marked A, B, C, D are all 2000 ticks and we calculated a 
potential for each range shown in A-1, B-1, C-1, D-1. In these potential figures, we set scales of 
horizontal axis and bottom axis the same size. We show the estimated potential curvature b for 
each potential. 
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FIG. 2: Reproduced potentials by deterministic dealer model. The ranges marked A', B', C', D' 
are 2000 ticks as well as Fig(TJWe change the foreseeing parameter d for each range. In range A', 
d is -0.5. In other ranges B', C, D' we change d=0, 0.5, 1.0 respectively. For other parameters we 
identically set as follows: N (the number of dealer) is 300, L (spread) is 1.0 and the range of q is 
[0.01,0.02]. 
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FIG. 3: The relation between potential curvature b* plotted horizontal axis and dealers' foreseeing 
effect d plotted bottom line. The b* means the average of tick time for b. In this simulation we set 
the parameters as follows: N (the number of dealer) is 300, L (spread) is 1.0 and the range of q 
is [0.01,0.02]. 
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